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were  fabricated  and  assembled  into  optical  processors;  and  the  perfor- 
mance of  the  processors  was  measured  and  compared  to  the  analysis.  Both 
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benefit  from  the  unique  properties  of  holographic  optics.  A useful 
quality  criterion  that  was  developed  for  analyzing  matched  filter  opti- 
cal processors  is  the  standard  deviation  phase  error  criterion, 
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Final  Scientific  Report 

OPTICAL  PROCESSORS  USING  HOLOGRAPHIC  OPTICAL  ELEMENTS 


This  report  describes  a two-phase  effort  to  investigate  holographic 
optical  elements  for  use  in  optical  processors.  In  the  first  phase, 
holographic  optics  for  a matched  filter  optical  processor  and  holographic 
optics  for  a Fourier  transform  optical  processor  were  designed  and  ana- 
lyzed. In  the  second  phase,  holographic  optics  were  fabricated,  tested, 
and  assembled  into  optical  processors;  optical  processing  experiments 
were  undertaken,  and  the  performance  of  the  processors  was  compared  with 
that  predicted  by  the  analysis.  The  results  of  this  research  suggest 
/ that  holographic  lenses  are  an  attractive  alternative  to  conventional 
optics  for  those  optical  processing  applications  that  benefit  from  the 
unique  advantages  of  holographic  optics.  These  advantages  include  their 
light  weight,  compact  design,  low  cost,  and  ability  to  perform  multiple 
functions. 

Major  results  of  this  research  are  described  in  Appendix  A of  this 
report,  which  is  prepared  in  the  form  of  a paper  entitled,  "Holographic 
Optics  for  a Matched-Filter  Optical  Processor,"  [1]  which  is  being  sub- 
mitted for  publication.  Another  publication  resulting  from  this  work  is 
Ref.  2. 

Important  topics  not  covered  in  Appendix  A comprise  the  remainder  of 
the  body  of  this  report.  The  topics  include  the  implementation  of  a real- 
time optical  processor  and  a s)Tithetic  aperture  radar  processor. 


1.  J.R.  Fienup  and  C.D.  Leonard,  "Holographic  Optics  for  a Matched-Filter 
Optical  Processor,"  submitted  for  publication  to  Applied  Optics. 

2.  J.R.  Fienup,  W.S.  Colburn,  B.J.  Chang,  and  C.D.  Leonard,  "Compact  Real- 
Time  Matched-Filter  Optical  Processor,"  Proc.  SPIE  118,  Optical  Signal 
and  Image  Processing,  p.  21  (August  1977). 
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REAL-TIME  OPTICAL  PROCESSOR 

The  optical  processor  used  for  the  experiments  described  in  Appen- 
dix A is  shown  in  Figure  1 (a  schematic  diagram  of  which  is  shown  in 
Figure  1(b)  of  Appendix  A).  A collimated  beam  illuminates  the  input 
transparency  (S) , the  Fourier  transform  of  which  is  brought  to  focus 
at  the  matched  filter  plane  (F)  by  the  reflective  holographic  lens  (H^). 
The  wavefront  emerging  from  the  filter  is  brought  to  focus  at  the  out- 
put plane  (OD)  by  a second  reflective  holographic  lens  (H2)-  The  colli- 
mating lens  L^  provides  a plane  wave  reference  beam  for  recording  the 
matched  filter o 

In  addition  to  those  described  in  Appendix  A,  correlation  experi- 
ments were  performed  with  a photoconductor-thermoplastic  device  as  the 
input  in  order  to  demonstrate  the  concept  of  a lightweight  and  compact 
real-time  matched-filter  optical  processor » The  photoconductor-thermo- 
plastic  device  was  developed  in  a parallel  effort  sponsored  by  the  Air 
Force  Office  of  Scientific  Research  under  Contract  NOo  F44620-76-C-0053. 
The  experimental  setup,  including  the  real-time  input,  is  shown  in  the 
diagram  of  Figure  2 and  the  photograph  in  Figure  3o  A transparency  (A) 
is  illuminated  with  incoherent  light  and  imaged  via  a beamsplitter  (B) 
onto  the  photoconductor-thermoplastic  device  (S)  which  is  in  the  input 
plane  to  the  optical  processor.  A grating  (G)  is  sandwiched  with  the 
transparency  (A)  in  order  to  provide  a spatial  frequency  offset  as  re- 
quired by  the  band-pass  response  of  the  photoconductor-thermoplat tic 
device.  The  photoconductor-thermoplastic  device  acts  as  an  incoherent- 
to-coherent  transducer  enabling  new  data  to  be  entered  in  real-time  at 
the  input  of  the  coherent  optical  processor.  Figure  4 shows  a scan 
through  correlation  output  achieved  with  the  real-time  processor.  A 
signal-to-noise  ratio  exceeding  30  dB  and  a 3 dB  spot  width  of  47  pm 
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Figure  2.  Real-time  matched-filter  optical  processor 
using  holographic  optics.  Dotted  lines 
indicate  imaging  with  incoherent  light, 
solid  lines  indicate  coherent  light  for 
optical  processing. 


I'iKuri'  3.  Real-tinif  maLi’lied-f  i 1 tor  optioal  prooossor  iisinp  lin  lo^raph  i o 

optioal  olomonts. 
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were  realized^  The  spot  size  is  very  close  to  the  theoretically  predicted 
diffraction-limited  size  of  44  ym  that  would  arise  from  the  7,1  mm  aperture 
used  in  the  filter  plane. 

In  this  case,  the  performance  of  the  holographic  optics  far  exceeded 
the  space-bandwidth  product  requirements  of  the  input  transducer,  which 
was  limited  to  200  x 200  picture  elements  due  to  the  frequency  of  the 
grating  sandwiched  with  the  transparency.  If  the  processor  were  made  to 
be  compact  (using  focal  lengths  of  100  mm  or  less) , then  the  performance 
of  the  holographic  optics  would  more  closely  match  that  of  the  input  trans- 
ducer, 

SAR  DATA  PROCESSOR 

Synthetic  aperture  radar  (SAR)  optical  data  processors  are  another 
potential  application  for  holographic  optics.  A holographic  lens  was  used 
for  Fourier  transform  processing  of  radar  film.  The  resulting  image  com- 
pares favorably  with  an  image  produced  by  a conventional  optical  processor. 

The  imagery,  which  is  classified  SECRET,  is  not  included  in  this  report. 

The  potential  light  weight  and  compact  design  of  holographic  optics  are 
attractive  features  for  optical  processors  for  airborne  radar  systems.  The 
predicted  performance  of  a simple  holographic  lens  for  two-dimensional  Fourier 
transformations  is  given  by  Figure  6 of  Appendix  A, 

CONCLUSIONS 

Both  analytical  and  experimental  results  obtained  under  this  effort 
show  that  holographic  optics  perform  well  as  the  Fourier  transforming 
lenses  in  optical  processors.  Images  of  moderate  space-bandwidth  product 


15 


can  be  processed  by  a very  compact  and  light  weight  optical  processing 
system  using  holographic  lenses  of  a very  simple  design^  If  combined 
with  a real-time  incoherent-to-coherent  transducer,  such  as  the  photo- 
conductor-thermoplastic device,  holographic  optics  can  be  used  to  form 
a terrain-matching  optical  correlator  suitable  for  terminal  guidance 
applicationSo  If  longer  focal  lengths  are  allowed,  then  holographic 
lenses  are  capable  of  processing  imagery  of  large  space-bandwidth 
product. 

The  results  of  the  analysis  of  the  holographic  lenses  performed  in 
this  study  was  generally  verified  by  the  experimental  results.  The 
standard  deviation  phase  error  criterion  developed  for  the  matched-filter 
optical  processor  was  found  to  be  a useful  guide  for  designing  and  analyzing 
holographic  Fourier  transforming  lenses. 

Further  improvement  in  the  performance  of  the  holographic  lenses 
could  come  in  three  areas:  the  use  of  multiple-element  Fourier  transform 
lenses,  the  use  of  holograms  and  output  detectors  on  curved  surfaces,  and 
the  use  of  holographic  lenses  recorded  with  nonspherical  wavefronts^ 
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Appendix  A 

HOLOGRAPHIC  OPTICS  FOR  A MATCHED-FILTER  OPTICAL  PROCESSOR 

by 

Jo  Ro  Fienup  and  Co  Do  Leonard 

Environmental  Research  Institute  of  Michigan 
P.  0.  Box  8618,  Ann  Arbor,  Michigan  48107 

ABSTRACT 

The  requirements  on  and  the  performance  of  holographic  Fourier 
transform  lenses  for  a matched-filter  optical  processor  are  discussedo 
The  holographic  aberrations  were  analyzed  using  a ray-tracing  computer 
program,  and  holographic  lenses  were  fabricated  and  assembled  into  an 
optical  processor.  Analysis  and  experimental  results  show  that  holo- 
graphic optics  are  capable  of  processing  images  of  large  space-bandwidth 
product  and  can  be  an  attractive  alternative  to  conventional  optics. 


(This  work  was  supported  by  the  Air  Force  Office  of  Scientific  Research 
under  Contract  No.  F44620-76-C-0047. ) 
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I.  INTRODUCTION 

Coherent  optical  processors  have  proven  useful  for  a number  of 
applications  [1],  The  unique  properties  of  holographic  optical  elements 
offer  a number  of  advantages  over  glass  lenses  for  use  in  optical  process- 
ors. Compared  with  conventional  optics,  holographic  lenses  are  thin  and 
light  weight.  Made  by  a photographic  process,  they  are  potentially  in- 
expensive. They  can  be  made  in  reflection  as  well  as  transmission,  allow- 
ing for  compact  folded  light  paths.  Additional  advantages  Include  the 
possibility  of  multiple  functions  in  the  same  area  of  a given  hologram 
(such  as  focusing,  beam  splitting,  and  spectral  filtering),  multiple  super- 
imposed holograms,  the  ability  of  a hologram  on  a curved  surface  to  function 
Independently  of  shape  of  the  surface,  and  the  ease  of  recording  holographic 
elements  with  non-spherical  wavefronts. 

Reflection  holograms  recorded  in  dichromated  gelatin  can  satisfy  the 
requirements  of  an  optical  processor  for  high  diffraction  efficiency,  high 
angular  bandwidth,  and  low  optical  noise  [2].  Dichromated  gelatin  holograms 
of  good  optical  quality  can  diffract  over  95%  of  the  incident  light  into  the 
first  diffracted  order  and  can  have  angular  bandwidths  of  several  degrees. 

We  analyzed  the  aberrations  of  holographic  optics  in  a matched  filter 
optical  processor,  and  fabricated  and  tested  an  optical  processor  using 
holographic  lenses.  We  found  holographic  lenses  to  perform  well,  and 
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conclude  that  they  are  an  attractive  alternative  to  conventional  lenses 
for  those  applications  that  can  benefit  from  their  unique  properties. 


II.  ANALYSIS  OF  OPTICS  FOR  A MATCHED-FILTER  PROCESSOR 

For  comparison,  a conventional  matched  filter  optical  processor 
and  its  holographic  optics  counterpart  are  shown  in  Figure  1(a)  and  (b), 
respectively.  As  previously  noted,  the  striking  properties  of  the  holo- 
graphic optical  processor  are  its  compact  folded  geometry  and  its  light 
weight.  Referring  to  Figure  1,  the  input  (signal)  transparency  S,  which 
has  amplitude  transmittance  Sj^(x),  is  illuminated  by  a plane  wave  from 
the  collimator  (or  L^) . We  use  one-dimensional  notation  for  simplicity; 
the  extension  to  two  dimensions  is  obvious.  In  the  ideal  case,  the  first 
Fourier  transforming  element  (or  L^)  causes  the  wavefront  S^(u)  = 9’{s,(x)}. 
the  Fourier  transform  of  s^(x),  to  appear  at  the  filter  plane.  The  matched 
filter  has  been  made  to  have  an  effective  transmittance  S*(u)  = (^ 
where  S2(x)  is  a reference  image  and  (*)  denotes  the  complex  conjugate. 

The  second  Fourier  transforming  element  (also  referred  to  as  the  re-imaging 
element)  (or  L2)  causes  the  Fourier  transform  of  S^(u)  S*(u),  the  wave- 
front  emerging  from  the  filter,  to  appear  at  the  output  plane  OD.  For 
simplicity,  assume  that  the  focal  lengths  of  the  two  Fourier  transforming 
elements  are  the  same.  The  amplitude  of  the  output  is  given  by 
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g(x)  = Sj^(x)  ★ S^Cx) 


s^(Xj^)  s*(x+x^)dx^ 


(1) 


which  is  the  desired  cross-correlation.  In  Eq.  (1),  constant  factors  that 
depend  on  optical  parameters  are  ignored.  If  the  input  image  is  simply  the 
reference  image  shifted  so  that  it  is  centered  at  x^,  then  the  output  be- 
comes the  autocorrelation 


g(x) 


s(x^)  s (x+x^+x^)dx^ 


(2) 


which  has  its  peak  at  -x  : 

o 


Ao  = g(-x^) 


(3) 


where  the  second  line  of  the  equation  above  results  from  Parseval's  tlieorem. 
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As  with  conventional  optics,  holographic  optics  have  aberrations 
that  limit  the  performance  of  optical  processors.  For  conventional 
optics,  aberrations  can  be  reduced  by  the  use  of  well-designed  multiple- 
element  lens  systems.  However,  a similar  reduction  of  aberrations  is 
usually  not  practical  with  holographic  optics.  Holographic  elements 
ordinarily  must  be  separated  by  a considerable  distance  in  order  to 
separate  the  desired  order  of  diffraction  from  the  spurious  ones, 

although  this  effect  is  minimized  by  volume  phase  holograms.  For  re- 
flection holograms,  the  element  following  the  hologram  cannot  be  allowed 

to  block  the  beam  incident  on  the  hologram.  Furthermore,  due  to  their 
off-axis  nature,  a group  of  hologram  lenses  would  be  more  difficult 

to  align  than  on-axis  optics.  Besides,  the  use  of  multiple-element 
holographic  optics  would  negate  the  advantages  of  light  weight  and  low 
cost.  For  these  reasons,  we  restricted  our  attention  to  simple  single- 
element holographic  lenses  on  plane  substrates. 

Unlike  a simple  lens  with  spherical  surfaces,  a holographic  lens 
recorded  with  spherical  wavefronts  will  introduce  no  phase  errors  if 
illuminated  with  the  wavefront  for  which  it  was  made.  However,  as  the 
curvature  and  field  angle  of  the  incident  wavefront  differ  from  the 

reference  beam  used  to  record  the  hologram,  aberrations  appear,  as 
analyzed  by  Champagne  [3].  The  chromatic  aberrations  are  not  a factor 
since  the  optical  processor  can  operate  with  monochromatic  light  of  the 
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same  wavelength  used  to  record  the  hologram.  Since  a hologram  lens  In 
an  optical  processor  operates  on  a complex  wavefront  comprised  of  the 
sum  of  many  wavefronts  of  various  field  angles,  it  necessarily  exhibits 
aberrations  [4].  The  nonchromatic  holographic  aberrations  are  the  limit- 
ing factor  in  the  performance  of  the  processor. 

A.  Standard  Deviation  Criterion 

In  order  to  specify  the  maximum  tolerable  aberrations  that  the 
Fourier  transforming  elements  can  be  allowed  to  introduce,  it  is  first 
necessary  to  determine,  in  general,  the  effects  of  aberrations  in  the 
Fourier  transforming  elements  on  the  correlation  output  of  the  matched 
filter  optical  processor.  However,  due  to  the  operation  of  the  matched 
filter  combined  with  the  space  variance  of  the  aberrations,  it  is  not 
possible  to  analyze  the  entire  system  by  conventional  wavefront  analysis 
or  ray-tracing  methods.  Therefore,  it  is  necessary  to  simplify  the 
analysis,  say,  by  treating  only  one  of  the  elements  at  a time.  In  order 
to  judge  a single  Fourier  transforming  element,  we  used  a standard  devi- 
ation error  criterion,  which  is  derived  below.  This  derivation  applies 
to  conventional  as  well  as  holographic  optics. 

Referring  again  to  Figure  1,  suppose  that  introduces  a phase 
error  <()^(x).  Strictly  speaking,  Is  a function  not  only  of  the  holo- 
gram coordinate,  but  also  of  the  coordinate,  u,  of  the  filter  plane; 
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that  Is,  is  space  variant.  We  will  continue  this  derivation  with  the 
simplit  ication  that  4)^^  is  space  invariant;  a justification  for  this  simpl 
fication  will  be  given  later.  A further  simplification  is  that  the  phase 
error  is  introduced  near  the  input  plane.  This  is  strictly  true  for 
certain  processor  geometries  and  approximately  true  for  slowly  varying 
phase  errors  for  the  geometry  of  Figure  1.  The  justification  of  the 
latter  simplification  relies  on  the  fact  that  slowly  varying  phase  wave- 
fronts  change  slowly  as  they  propagate.  Since  the  phase  errors  involved 
are,  in  fact,  small  and  slowly  varying,  this  is  a good  approximation. 

Thus,  for  a slowly  varying  phase  error  (j)^,  the  wavefront  arriving  at  the 
filter  is  proportional  to 


[S(u)exp(-j  ux^)]  * ^ I exp(j(})^)(- 


where  X is  the  wavelength  of  light,  f is  the  focal  length  of  the  Fourier 
transforming  element,  is  the  translation  of  the  center  of  the  input 

from  the  optical  axis,  and  ¥:  denotes  convolution.  This  wavefront,  multi 
plied  by  the  filter  function  S*(u),  is  approximately  an  apodized  plane 
wave  (if  no  phase  errors  were  present,  its  phase  would  be  -2ti  ux^/Xf)  . 

It  is  Fourier  transformed  and  brought  to  focus  in  the  output  plane 
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by  which  introduces  a slowly  varying  phase  error  <t)2  (which  depends 
on  Xq).  The  amplitude  of  the  wavefront  at  the  output  is  given  by 


g(x)  = ([s(u)exp(-_i  Iy  ux^)]*  ^ j exp(j  ) S*  (u)exp(j  t})^)  | (5) 


where  aberrations  introduced  by  the  matched  filter  are  ignored  [5].  if 
no  phase  errors  are  present,  <J)^  = ())^  = 0,  then  Eq,  (5)  would  reduce  to 
Eq.  (2). 

The  effects  of  phase  errors  (p^  are  made  more  clear  by  setting  4)^^  = 0 
in  Eq.  (5),  which  then  reduces  to 


g(x) 


O 

I 


2 Pit 

[S(u)|  exp[j  ^^(u) ]exp[-j  ^u(x+x^)]  du  (6) 


Evaluating  at  the  central  peak,  x = the  ratio  of  the  intensity  of 

the  central  peak  to  its  theoretical  maximum  value  is 
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where  is  a weighted  average  of  the  k power  of  (p^(u).  The 

2 

j weighting  function  ls(u)|  is  simply  the  intensity  pattern  of  the  wave- 

I 


front  emerging  from  the  spatial  filter.  For  reasons  described  later. 
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out  of  the  system,  linear  phase  errors  should  not  be  included  in  the 
calculation  of  O^. 

The  effects  of  phase  error  (J)^,  introduced  by  the  first 
Fourier  transforming  element  are  made  more  clear  by  setting  0^=0 

in  Eq.  (5),  which  then  reduces  to 


g(x) 


oo 


s (C)  exp  [ j<})^(C+x^)  ] s*  (5+x^-Hc)d5 


(10) 


Evaluating  at  the  peak,  x = “X^,  this  becomes 


g(-Xo)  = 


s(C) I ^exp[j$(C+x^) ]dC 


(11) 


Again  expanding 
weighted  sum  of 


the  exponential  and  interpreting  the  integral  as  a 
the  phase  errors,  we  find  that  for  small  phase  errors. 


R == 


(12) 
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The  approximations  assumed  in  the  derivation  of  Eqo  (5)  were  that 
be  small,  slowly  varying,  and  space-invar iantc  While  generally 
will  be  small  and  slowly  varying,  it  will  not  be  space-invariant , In 
most  cases,  increases  with  increasing  field  angles,  or,  equivalently, 
with  increasing  |u|  in  the  filter  plane,  and  is  maximum  for  u coordinates 
at  the  edge  of  the  filter  plane,.  The  justification  for  the  space-invariance 
approximation  is  as  follows.  Since  the  purpose  of  the  matched  filter  optical 


processor  is  to  produce  a sharp  peak  in  the  output  plane,  it  is  the  value  of 
for  points  near  the  edge  of  the  filter  plane  that  determines  the  perform- 
ance of  the  processor.  Using  0^  at  the  edge  of  the  filter  is  appropriate, 
since  the  information  at  the  edges  of  the  matched  filter  is  ordinarily  empha- 
sized in  order  to  increase  the  selectivity  of  the  filter  and  increase  the 
sharpness  of  the  autocorrelation  peak  [6],  This  is  usually  unavoidable 


t 
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since  it  is  difficult  to  record  the  large  dynamic  range  of  S*(u)  in  the 
matched  filter;  but  it  is  also  highly  desirable,,  The  wavefront  emerging 
from  the  matched  filter  has  an  intensity  distribution  that  is  more  nearly 
constant  or  of  an  annular  shape,  rather  than  being  sharply  peaked  in  the 
center  as  |S(u)|  would  beg  In  this  case,  the  output  is  no  longer  truly 
an  autocorrelation;  it  more  closely  resembles  a point  imageo  Thus,  we 
are  led  to  base  our  analysis  of  the  performance  of  element  on  the  value 
of  for  field  angles  corresponding  to  the  edges  of  the  filter^ 


I 


t 


A further  justification  of  basing  the  analysis  on  at  the  edge 
of  the  filter  is  as  followso  Suppose  that  for  a given  input  size  and 
system  design,  goes  beyond  a certain  critical  level  (the  value  of 
that  critical  level  will  be  shown  later)  for  positions  in  the  filter 
plane  with  |u(>  Then  an  aperture  should  be  placed  at  the 

filter  plane  to  block  out  the  portion  of  the  wavefront  for  which  |u|> 
|Umajji»  since  that  portion  of  the  wavefront  will  lead  to  a degradation 
rather  than  a reinforcement  of  the  autocorrelation  peak  at  the  processor 
outputs  That  value  of  |u  I defines  the  size  of  the  filter  used  and, 
consequently,  the  maximum  spatial  frequency  component  of  the  input  trans- 
parency that  is  usedg  Thus,  the  maximum  usable  space-bandwidth  product 
of  the  input  and  the  resolution  at  the  output  is  determined  by  the  value 
of  |u  I for  which  0,  goes  above  the  critical  level » One  can  therefore 
ignore  the  space  variance  of  the  phase  error  4)^^  and  concentrate  on  the 
value  of  as  computed  for  the  edge  of  the  matched  filters 
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To  test  the  validity  of  Eq.  (12)  and  to  more  accurately  establish 

the  effects  of  phase  errors  we  digitally  computed  the  aberrated 

autocorrelation  integral  of  Eq.  (5)  with  = 0 for  the  hypothetical 

one-dimensional  object  shown  in  Figure  2(a)  with  various  phase 

errors  (f)^.  Figure  2(b)  shows  the  central  portion  of  the  autocorrelation 

output  with  no  phase  error  (uppermost  curve)  and  with  = 1/16  wavelength 

n*'^-order  phase  error  for  n = 1 to  4,  respectively.  That  is,  one  curve 

results  from  a linear  phase  error  of  .2165  wavelengths  peak-to-peak,  the 

second  from  a quadratic  phase  error  of  .2097  wavelengths,  etc.  The  peaks 

of  all  four  curves  fall  into  the  range  of  0.842  to  0.868,  in  very  good 

2 

agreement  with  the  value  1 - (27r/16)  = 0.846  predicted  by  Eq.  (12). 

Figure  2(c)  shows  the  same  thing  with  = 1/8  wavelength.  In  this  case, 
the  peaks  fall  in  the  range  0<,478  to  0.567,  not  in  very  good  agreement 
with  the  value  0.383  predicted  by  Eq.  (12).  Figure  2(d)  shows  the  same 
thing  with  = 1/4  wavelength.  From  Figures  2(c)  and  2(d),  it  is  apparent 
that  Eq.  (12)  is  no  longer  reliable  for  > 1/8  wavelength.  However,  for 
much  greater  than  1/8  wavelength  the  autocorrelation  peak  would  no  longer 
be  detectable. 

From  the  above  results,  we  can  formulate  the  following  standard  devi- 
ation phase  error  criterion:  the  standard  deviation  of  the  total  phase 
error  should  be  kept  to  1/8  wavelength  in  order  to  insure  that  the  auto- 
correlation peak  remain  easily  detectable.  This  corresponds  to  the  follow- 
ing peak-to-peak  errors  (in  wavelengths):  .433  of  linear,  .419  of  quadratic. 
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.660  of  cubic,  and  ,468  of  fourth-order.  As  shown  in  Figure  2(c),  = 

1/8  wavelengths  results  in  a greatly  reduced,  but  still  detectable  , 
autocorrelation  peak;  and  as  shown  in  Figure  2(d),  = 1/4  wavelengths 

results  in  an  autocorrelation  peak  that  would  be  lost  in  the  noise. 

As  discussed  earlier,  it  is  desirable  to  record  the  matched  filter 
in  a nonlinear  way  so  as  to  block  the  lower  and  emphasize  the  higher 
spatial  frequencies.  In  a computer  simulation,  this  effect  was  approxi- 
mated in  Eq.  (5)  by  multiplying  S*(u)  by  a high-pass  stop  of  the  form 
[1-rect (u/u^) ] , As  expected,  this  resulted  in  an  output  peak  that  is 
much  sharper  and  has  a higher  signal-to-noise  ratio.  An  example  of 
this  is  shown  in  Figure  2(e),  for  which  = 0 and  the  high-pass  stop 

was  20%  of  the  total  filter-plane  aperture  width.  Figure  2(f)  shows  the 

til 

results  with  = 1/4  wavelength  n order  phase  error  for  n = 1 to  4, 
respectively.  Since  the  output  peaks  are  easily  detectable  above  the 
noise,  these  curves  suggest  that  the  emphasis  of  the  higher  spatial 
frequencies  also  serves  to  make  the  process  less  sensitive  to  the  slowly 
varying  phase  errors.  Thus,  the  standard  deviation  phase  error  quality 
crieterion  developed  above  for  the  case  of  linear  filter  recording  should 
be  taken  as  a lower  bound  on  the  performance  of  the  Fourier  transforming  lens. 

B.  Aberration  Computation 

Champagne's  equations  [3]  can  be  used  to  compute  the  third-order 
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\ holographic  aberrations,  but  they  are  inconvenient  for  computing  the 


standard  deviation  phase  error,  O,  Furthermore,  they  do  not  include 
higher-order  aberrations  without  adding  considerable  complexity.  Con- 
sequently, to  accurately  and  efficiently  analyze  and  optimize  the  de- 
sign of  the  holographic  optics,  we  used  a holographic  ray-tracing  com- 
puter program  called  the  Holographic  Optics  Analysis  and  Design  (HOAD) 
program,  developed  at  ERIM  [7], 

The  hologram  H shown  in  Figure  1(b)  is  a collimator  and  is  re- 
c 

constructed  in  exactly  the  same  way  that  it  was  recorded.  Consequently, 
it  has  no  holographic  aberrations  and  produces  an  ideal  plane  wave. 

In  analyzing  H^,  the  first  Fourier  transforming  element,  it  is 
important  to  note  that,  whatever  aberrations  H^  introduces,  if  the  same 
element  is  used  to  record  the  matched  filter,  then  it  is  a matched  filter 
of  an  aberrated  reference  image.  If  the  untranslated  reference  image  is 
used  as  the  input,  and  the  same  aberrations  are  introduced,  then  it  should 
correlate  perfectly.  This  would  seem  to  indicate  that  a very  poor  lens 
could  be  used  for  H^^,  and  indeed,  that  is  the  case  if  the  image  is  not 
translated  in  the  input.  However,  as  the  input  is  translated,  the  effect- 
ive aberrations  change,  and  the  correlation  peak  degenerates.  Thus,  the 
phase  error  of  interest  is  not  the  absolute  phase  error,  but  is  the  dif- 
ference between  the  phase  error  when  the  center  of  the  input  image  is 
translated  to  x^  and  the  phase  error  when  the  center  of  the  image  is  on 
the  optical  axis  (i.e.,  as  the  matched  filter  was  recorded).  Consequently, 
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some  of  the  requirements  for  the  Fourier  transforming  lens,  for  example, 
that  the  Abbe  sine  condition  must  be  fulfilled  [8],  may  be  relaxed. 
Generally,  the  phase  error  differences  are  minimized  when  the  absolute 
phase  error  is  minimized,  and  low-aberration  lenses  are  required  for  the 
matched  filter  optical  processor. 

To  analyze  H^,  we  considered  one  plane-wave  spatial  frequency  com- 
ponent of  the  input  at  a time,  by  tracing  crossed  fans  of  parallel  rays 
at  a given  field  angle  from  the  input  plane,  through  H^,  to  the  filter 
plane.  Often  the  translation,  x^,  of  the  input  image  relative  to  the 
reference  image  is  comparable  in  size  to,  or  larger  than,  either  the 
input  image  or  the  reference  image  (in  many  applications,  one  of  the 
two  images  covers  a much  smaller  area  than  the  other  — we  refer  to  the 
smaller  of  the  two  as  the  "image  of  interest").  In  that  case,  the  pre- 
dominant term  in  the  difference  of  the  phase  errors  is  a linear  term. 

That  is,  the  predominant  effect  is  that  if  a plane-wave  component  of  the 
image  comes  to  focus  at  u^  when  the  image  is  centered  on  the  optical  axis, 

then  it  will  come  to  focus  at  u,  + Au  when  the  image  is  centered  at  x ; 

1 o 

and  the  linear  phase  term  is  proportional  to  Au.  The  difference  of  the 
phase  errors  can  be  approximated  by  computing  the  aberrations  of  the  plane- 
wave  component  of  the  input  image  with  respect  to  the  position  u^^  in  the 
filter  plane  (rather  than  computing  the  aberrations  with  respect  to  its 
own  image  point  at  u^  + Au) . 
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The  analysis  of  Is  simpler,  since  it  only  brings  a quasi-plane 
wave  to  focus  at  the  output,  and  linear  phase  errors  are  not  considered 
significant.  Because  of  this,  the  requirements  on  are  not  as  strin- 
gent as  on  H^,  nor  are  they  as  stringent  as  they  would  be  for  some  other 
types  of  processors. 


C.  Optimization  and  Analysis 


Holographic  optics  has  a somewhat  different  set  of  parameters  for 
design  optimization  than  is  available  to  conventional  optics.  The  over- 
all geometry  of  a holographic  optical  system  can  be  adjusted  by  appropriate 
choices  of  focal  lengths,  offset  angles,  and  transmissive  vs.  reflective 
holograms.  Since  the  gross  imaging  properties  of  an  individual  holographic 
element  do  not  depend  on  the  angular  orientation  or  shape  of  its  surface, 
the  angular  orientation  of  the  element  with  respect  to  the  rest  of  the 
system  and  the  curvature  of  the  surface  of  the  hologram  substrate  can  be 
designed  to  minimize  aberrations.  A hologram  can  be  made  at  one  wavelength 
and  read  out  at  a second  wavelength.  The  thickness  and  refractive  index 
modulation  of  a holographic  element  can  be  manipulated  to  produce  a desired 
Bragg  angular  and  wavelength  selectivity.  Analogous  to  conventional  lenses, 
holographic  elements  of  a given  focal  length  can  have  different  shape  factors, 
depending  on  the  curvatures  of  the  reference  and  object  beams  used  to  record 

the  hologram. 

In  general,  aberrations  tend  to  be  minimized  when  the  offset  angles 
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are  minimized,  so  the  minimum  offsets  allowed  by  the  system  geometry  were 
used.  It  was  assumed  that  the  processor  operates  at  the  same  wavelength, 

514.5  nm,  at  which  the  holograms  are  recorded,  which  also  tends  to  mini- 
mize aberrations.  In  order  to  keep  the  system  as  simple  and  economical 
as  possible,  only  single  elements  on  plane  surfaces  were  considered. 

1.  Before-the-Lens  Processor 

Using  the  HOAD  program,  we  optimized  the  shape  factors  of  and 
H2  and  their  angular  orientations  with  respect  to  the  rest  of  the  sys- 
tem. As  depicted  in  Figure  1(b),  and  the  filter  should  be  parallel 
to  one  another  and  normal  to  the  axis  between  them.  Similarly,  H2  and 
the  output  plane  should  be  parallel  to  one  another,  and  normal  to  the 
axis  between  them.  (Note  that  the  input  plane  S is  parallel  to  H^,  as 
shown,  making  it  tilted  with  respect  to  the  S-H^^  axis.  If  S is  rotated 
to  be  normal  to  the  S-H^  axis,  then  the  output  plane  OD  rotates  to  be- 
come parallel  to  S,  but  tilted  with  respect  to  the  H2-OD  axis.)  The 
recording  geometry  for  is  shown  in  Figure  3 for  a 10°  offset,  where 
the  point  source  for  the  object  beam  is  located  at  the  center  of  the 
filter  plane  and  the  reference  beam  is  a plane  wave  conjugate  to  the 
S-Hj^  axis.  The  recording  geometry  of  is  the  mirror  image  of  Figure  3, 
where  the  point  source  for  the  object  beam  is  located  at  the  center  of 
the  output  plane  and  the  reference  beam  is  a plane  wave  conjugate  to  the 
filter  — H2  axis. 
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The  HOAD  program  was  used  to  compute  a^,  the  standard  deviation 
phase  error  of  H^,  as  a function  of  field  angle  (or  filter  plane  lo- 
cation), diameter  D of  the  image  of  interest,  translation  x of  the 
s o 

input,  and  focal  length.  The  performance  can  be  determined,  for  ex- 
ample, by  fixing  all  the  parameters  except  x^,  then  increasing  x^  until 
equals  1/8  wavlength,.  This,  along  with  the  diameter  of  the  image  of 
interest,  defines  D = + 2x^,  the  greatest  width  of  an  image  that  can 

be  processed,  that  is,  within  which  the  image  of  interest  can  be  located 
and  still  produce  a detectable  correlation  peak.  Combined  with  the  field 
angle  (corresponding  to  the  maximum  spatial  frequency  component  of  the 
image) , the  width  gives  the  maximum  space-bandwidth  product  (number  of 
picture  elements)  that  can  be  processed  without  a loss  of  resolution. 

Two  families  of  operating  curves  for  a 100mm  focal-length  are 

shown  in  Figure  4.  The  varying  parameters  for  the  different  curves  are 

D and  the  maximum  spatial  frequency,  respectively.  These  curves  may  be 
s 

used  as  follows.  If  it  is  required  that  the  resolution  of  the  matched 
filtering  operation  be  equivalent  to  20  cyc/mm  in  the  input  plane,  and 
that  the  image  of  interest  be  detectable  anywhere  within  a total  image 
area  of  750  x 750  picture  elements,  then  the  image  of  Interest  could 
contain  no  more  than  100  x loO  picture  elements.  Two  facts  made  evident 
by  Figure  4 are  that  the  performance  improves  if  ttie  image  is  demagnified 
(i.e.,  it  performs  better  for  higher  spatial  frequencies,  but  proportionally 
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smaller  image  widths),  and  that  there  is  a trade-off  between  and  D 
(i.e.,  the  larger  the  image  of  interest,  the  less  it  can  be  translated). 

The  aberrations  of  the  second  Fourier  transforming  element,  H^, 
were  also  computed.  Since  it  was  found  that  the  aberrations  of 
generally  were  the  limiting  factor  in  the  performance  of  the  system, 
the  less  interesting  results  on  H2  will  not  be  given  here. 

A most  important  parameter  is  the  focal  length,  F,  of  the  Fourier 

transforming  elements.  Using  Champagne's  equations  [3],  it  can  be  shown 

that,  for  a given  image  size  and  field  angle,  third-order  coma  is  pro- 
2 

portional  to  1/F  , and  third-order  astigmatism  is  proportional  to  1/F. 

Consequently,  the  performance  of  the  processor  improves  greatly  with 

increasing  focal  lengths.  As  an  example,  Figure  5 shows  the  translation 

that  results  in  = 1/8  wavelength  vs.  focal  length,  for  various 

values  of  D . A highest  spatial  frequency  of  40  cyc/mm  and  an  offset  of 
s 

10°  was  assumed.  In  this  case,  10mm  is  equivalent  to  800  picture  ele- 
ments. Figure  5 indicates  that  a compact  processor  using  elements  with 
short  (100mm)  focal  lengths  can  process  imagery  of  a few  hundred  pic- 
ture elements  on  a side;  and  if  long  (>  1 meter)  focal  lengths  are 
allowed,  then  the  holographic  optics  can  handle  imagery  having  thousands 
of  picture  elements  on  a side. 

2.  Af ter-the-Lens  Processor 

An  alternative  to  the  configurations  shown  in  Figure  1 is  to  locate 
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the  input  plane  in  a converging  beam  after  a transmissive  first  Fourier 
transforming  lens  [5] , which  would  allow  for  easy  changes  in  scale  by 
translating  the  input  along  the  optical  axis.  In  this  case,  the  first 
Fourier  transforming  element,  which  could  be  combined  with  the  colli- 
mator, would  produce  a perfect  converging  spherical  wavefront.  Never- 
theless, the  resulting  Fourier  transform  would  be  aberrated  (in  addition 
to  having  a spherical  phase  factor  associated  with  the  Fourier  transform). 


In  order  to  analyse  the  aberrations  of  the  af ter-the-lens  Fourier 

transform  geometry,  consider  a single  spatial  frequency  component  of  the 

input.  For  this  analysis,  we  use  Champagne's  notation  and  equations  [3], 

in  which  R is  the  radial  distance  of  a recording  or  reconstruction  point 

source  from  the  hologram,  a and  3 are  field  angles,  and  subscripts  R,  0, 

C,  and  I refer  to  reference,  object,  reconstruction,  and  image  beams, 

respectively.  The  input  then  can  be  modeled  as  a constant-frequency 

holographic  grating  with  recording  parameters  R = R = + °°,  a = 3„=  0, 

K 0 R R 

and  Qq  and  3q  are  the  angles  corresponding  to  the  spatial  frequency  under 
consideration.  For  simplicity,  consider  3q  = 0,  since  the  performance  will 
be  s}Tnmetrlc  in  the  x and  y dimensions.  The  reconstruction  beam  is  a spheri- 
cal wave  converging  to  the  center  of  the  matched  filter  plane,  and  is  given 
by  R = -F  and  a = 3^  = 0*  Then  Champagne's  equations  for  the  image 

c c ^ 
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location  reduce  to 


Rj  - Rc  - -f 


sin  = sin  a^,  or,  = a 


cos  a^.  sin  = cos  sin  = 0 ; 


and  the  coefficients  of  the  third-order  aberrations  reduce  to 


(13) 


(14) 


(15) 


"c  h 


sin  a. 


sin  a 


(16) 


(17) 
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and 


A = A =0  (20) 

y xy 


where  the  coefficients  S,  C,  and  A are  for  aberrations  of  the  spherical, 
comatic,  and  astigmatic  types,  respectively. 

It  is  interesting  to  compare  these  aberrations  with  those  of  the 

before-the-lens  processor.  For  the  bef ore-the-lens  processor,  using  the 

holographic  element  shown  in  Figure  2,  the  recording  parameters  are  R = 

R_  = F,  a„  = 10°,  and  6 = a-  = 3^  = 0;  and  for  the  reconstruction,  R = °°, 

UR  KUO  C 

and  and  3q  are  variable.  If  3q  is  set  to  zero  and  sin  a is  substituted 

for  sin  ol„  - sin  a , then  it  can  be  shown  that  the  aberration  coefficients 
0 R 

of  the  optimized  holographic  before-the  lens  processor  are  exactly  the  same 
as  those  given  above  for  the  af ter-the-lens  processor.  The  aberrations  in 


f 
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the  two  cases  are  not  exactly  the  same,  however,  since  a given  plane  wave 
spatial  frequency  component  of  the  input  to  the  before-the-lens  processor 
arrives  at  off  center  by  a distance  F sin  a„  Nevertheless,  using  the 
HOAD  program  it  was  found  that,  for  the  range  of  parameters  of  interest, 
the  performance  of  the  after-the-lens  processor  is  nearly  identical  to 
that  of  the  holographic  before-the-lens  processor „ 

In  addition  to  allowing  scale  changes,  the  after-the-lens  processor 
allows  for  a distance  of  F between  the  input  and  filter  planes,  as  com- 
pared to  2F  for  the  before-the-lens  processor  (this  advantage  is  diluted, 
however,  since  the  before-the-lens  processor  can  be  folded  to  reduce  the 
2F  length  back  to  F) „ Furthermore,  for  longer  focal  length  systems  for 
which  aF  is  comparable  to  D,  the  Fourier  transforming  element  in  the 
after-the-lens  processor  can  be  smaller  than  that  of  the  before-the-lens 
processor o A possible  disadvantage  of  the  after-the-lens  processor  is 
that  since  the  matched  filter  is  required  to  record  an  additional  spheri- 
cal phase  factor,  the  matched  filter  Itself  will  introduce  aberrations 
not  present  in  the  matched  filter  made  for  a before-the-lens  processor  [5,9] o 

3o  Fourier  Spectrum  Analysis 

The  holographic  Fourier  transforming  lenses  were  also  analyzed  for 
use  in  simple  Fourier  transform  optical  processors,  for  performing  Fourier 
spectrum  analysis  of  an  input  transparency.  Such  a processor  would  be  similar  to 
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the  first  half  of  the  matched-filter  optical  processor,  with  the  output 
plane  being  located  in  place  of  the  matched  filter » For  this  analysis, 
the  distortion  of  the  output  (linear  phase  errors)  was  ignored,  and  it 
was  assumed  that  the  optimized  holographic  lens  shown  in  Figures  3 with 
a 10°  offset  angle  is  used,  and  it  was  required  that  a not  exceed  1/8 
wavelength  at  40  cyc/mm  in  the  input.  Figure  6 shows  the  predicted 
linear  space-bandwidth  product  of  the  input  that  could  be  Fourier  trans- 
formed (which  is  proportional  to  the  maximum  allowed  width  of  the  input) 
as  a function  of  the  focal  length  of  the  holographic  lens. 


III.  HOLOGRAM  FABRICATION 

Dichromated  gelatin  plates  were  prepared  by  fixing  the  silver 
halide  out  of  5 x 5 inch  Kodak  649-0  microflat  plates,  then  sensitizing 
with  ammonium-dichromate  solution.  The  processing  procedures  used  were 
those  of  Chang  and  Colburn  [2],  which  are  similar  to  those  of  Chang  [10]. 
Due  to  the  physical  sizes  of  the  plate  holders  and  other  available  equip- 
ment, we  chose  the  focal  lengths  and  offset  angles  of  the  two  Fourier 

transforming  lenses  to  be  500  mm  and  20°,  respectively.  Using  the  514.5  nm 

2 

wavelength  from  an  argon-ion  laser,  typical  exposures  were  at  1 mW/cm  for 
eight  minutes. 

The  emulsion  side  of  the  plate  was  placed  away  from  the  point  source 
during  recording  and  away  from  the  input  and  filter  planes  when  in  the 
processor.  Consequently,  the  input  wavefront  passes  through  the  glass 
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substrate  both  before  and  after  reflection  from  the  holographic  lens» 

The  spherical  aberration  introduced  by  the  glass  substrate  into  the 
diverging  recording  beam  is  cancelled  by  the  opposite  aberration  intro- 
duced by  the  substrate  into  the  converging  reflected  wavefront  upon 
readouts  An  advantage  of  this  orientation  of  the  emulsion  is  that  the 
reconstructed  beam  is  reflected  without  passing  through  the  outer  sur- 
face of  the  emulsiono  Then  the  dichromated  gelatin  hologram  can  be 
sealed  (to  protect  the  emulsion  primarily  from  the  effects  of  humidity) 
by  cementing  on  a thin  glass  cover  plate,  without  having  to  consider 
the  optical  quality  of  the  cover  plateo 

The  diffraction  efficiency  of  a typical  holographic  lens  made  for 
the  experiments  is  70%„  Heavier  exposures  resulting  in  higher  diffraction 
efficiencies  were  avoided  since  they  also  result  in  the  recording  of  spur- 
ious holograms^  The  spurious  holograms  arise  from  beams  internally  re- 
flected from  the  air-glass  and  air-gelatin  surfaceso  For  the  highest 
quality  results,  those  unwanted  reflections  can  be  eliminated  by  record- 
ing the  hologram  in  a liquid  gate  with  antireflection  windows o 

The  diffraction  efficiency  of  the  holographic  lenses  must  remain 
high  for  a range  of  angles  of  illuminationo  For  the  first  Fourier  trans- 
forming lens,  the  required  angular  bandwidth  is  + a^,  where  is  the 
angle  corresponding  to  the  maximum  spatial  frequency  of  the  input;  for 
example,  ct  = lol8°  for  a frequency  of  40  cyc/mm„  For  the  second  Fourier 
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transforming  lens,  the  required  angular  bandwidth  is  + arctan  (x^/F), 
where  is  the  largest  allowed  translation  of  the  input;  for  example, 
if  Xjij/F  = 25 mm/500  mm,  then  the  required  angular  bandwidth  is  + 2o86°o 
The  angular  bandwidths  of  the  holographic  lenses  were  made  to  easily 
satisfy  these  requirements,  as  can  be  seen  from  the  curves  shown  in 
Figure  7 of  diffraction  efficiency  versus  angle  of  illumination  for  the 
central  portions  of  two  different  holographic  lenses.  Both  holographic 
lenses  were  recorded  with  20°  offset  angles  and  are  designed  to  operate 
at  angles  about  20°.  The  double  peaks  in  curve  A arise  from  the  symmetry 
of  the  Bragg  condition  about  the  angle  of  approximately  10°  to  which  the 
fringe  planes  are  normal.  The  right  peak  of  curve  A is  at  an  angle  greater 
than  20°  due  to  an  increase  in  the  thickness  of  the  hologram  emulsion.  A 
decrease  in  the  emulsion  thickness  caused  the  peaks  of  curve  B to  move 
together  to  form  a single  broad  peak.  The  emulsion  thickness  is  controlled 
primarily  by  the  concentration  of  the  sensitizer.  The  thickness  can  then 
be  fine-tuned  during  the  drying  process  in  which  the  holograms  are  baked 
in  an  oven  under  vacuum,  which  drives  moisture  from  the  emulsion.  The 
emulsion  thickness  can  be  decreased  by  additional  baking.  If  the  emulsion 
is  baked  excessively  and  becomes  too  thin,  then  it  can  be  made  thicker 
by  immersing  it  in  water  and  baking  it  again  to  a lesser  degree . In  this 
way  the  peak  of  the  angular  selectivity  curve  can  be  adjusted  as  desired , 
within  certain  limits. 

r 
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The  high  angular  bandwldths  of  the  holographic  lenses  result  from 
large  refractive  index  modulations,  determined  to  be  approximately  O.Olo 
Higher  modulation  could  be  achieved,  but  at  the  expense  of  optical  quality. 

Our  primary  consideration  was  the  phase  aberrations  of  the  holographic 
lenses o Before  preparing  them  for  exposure,  the  best  samples  from  a few 
boxes  of  649-0  plates  were  selected,  based  on  examination  of  flatness 
using  a Mach-Zehnder  interferometer o Since  washing  and  sensitization  of 
the  plates  was  done  after  the  interferometric  test,  it  did  not  matter  that 
the  plates  were  exposed  to  light  during  the  test„  After  being  recorded 
and  processed,  the  holograms  were  tested  for  aberrations  using  the  inter- 
ferometer geometry  shown  in  Figure  8o  The  position  of  the  point  source 
illuminating  the  hologram  was  the  same  as  that  used  to  record  it  and  so 

there  are  no  theoretical  holographic  aberrations.  The  aberrations  that 
appear  in  this  configuration  (the  "non-holographic  aberrations")  are  due, 
primarily,  to  variations  in  emulsion  thickness  and  edge  effects.  Figure  9 
shows  interf erograms  of  diameter  10  cm  of  two  different  holographic  lenses 

used  in  our  experiments . The  peak-to-peak  phase  error  of  lens  A is  about 
A/4  over  a 4 cm  diameter  and  A/2  over  an  8 cm  diameter,  and  of  lens  B is 
about  A/4  over  a 4 cm  diameter  and  about  A/2  over  a 6 cm  diameter.  These 
aberrations  are  sufficiently  low  for  the  purposes  of  our  experiments. 

IV.  CORRELATION  EXPERIMENTS 

The  holographic  lenses  described  in  the  previous  section  were 


A-29 




FORMEMLV  WILLOW  P»L>N  LABORATORIES.  THE  UNIVERSITY  OF  MICHIGAN 


assembled  into  a matched  filter  optical  processor,  and  matched  filter- 
ing experiments  were  performed^  Matched  filters  on  4 ^ 5 inch  Kodak 
131-01  microflat  plates  were  recorded  with  the  image  centered  in  the 
input  plane  and  using  to  perform  the  Fourier  transform;  referring 
to  Figure  1(b),  the  reference  beam  was  a plane  wave  along  the  filter- 
to-H2  axiso  The  intensities  of  the  reference  and  signal  beams  were 
adjusted  so  that  the  fringe  contrast  was  highest  for  spatial  frequencies 
above  15  cyc/mmo  This  resulted  in  matched  filters  that  are  opaque  in 
their  central  regions  because  of  heavy  overexposure.  The  filters,  giving 
heavy  weighting  to  the  higher  spatial  frequencies,  resulted  in  very 
sharply  peaked  "correlation"  outputs  [6]. 

In  one  experiment,  the  matched  filter  was  recorded  of  a 25mm 
diameter  area  of  the  input  image  shown  in  Figure  10.  The  filter  was 
limited  to  a 13.5mm  diameter  aperture,  corresponding  to  a field  angle 
of  0.77°  and  a spatial  frequency  of  26  cyc/mm.  The  central  3.5mm 
diameter  area  of  the  filter  was  opaque,  resulting  in  an  active  area  of 
the  filter  having  an  annular  shape.  A photometric  scan  through  the 
correlation  peak  for  x^  = 4 mm  is  shown  in  Figure  11.  The  half-power 
width  of  the  peak  shown  is  about  31  ym  and  it  is  more  than  30  dB  above 
the  noise.  The  sldelobe  ring  around  the  main  peak  is  a characteristic 
of  the  annular  shape  of  the  matched  filter.  The  intensity  of  the 
correlation  peak  as  a function  of  input  translation  x^  is  shown  in 
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Figure  12.  Even  at  - 21  ram,  for  which  its  intensity  is  15%  of  its 
value  for  = 0,  the  peak,  though  aberrated,  is  easily  detectable  and 
is  28  dB  above  the  noise,  as  shown  in  Figure  13. 

The  perforraance  of  the  holographic  optics  indicated  by  the  results 
above  is  somewhat  better  than  that  predicted  by  the  standard  deviation 
phase  error  criterion^  For  an  input  of  diameter  25mm  and  a field  angle 
of  0.- 77°,  the  HOAD  program  predicts  = 1/8  wavelength  for  x^  = 10  mm; 
but  the  measured  peak  Intensity  does  not  drop  to  50%  until  x^  = 13.5  mm, 
and  is  easily  detectable  for  x^  beyond  21  mm.  This  discrepancy  is  due  to 
the  inaccuracy  of  the  space  invariance  assumption  and  the  decreased  sensi- 
tivity of  the  processor  to  slowly  varying  phase  errors  when  the  higher  spa- 
tial frequency  information  is  emphasized.  Thus,  the  standard  deviation  phase 
error  criterion  should  be  taken  as  a lower  bound  on  the  performance  of  the  processor. 

V.  IMAGING  EXPERIMENT 

Another  test  of  the  holographic  Fourier  transform  lenses  was  to  use 
one  as  the  first  element  in  a two-lens  imaging  system,  as  shown  in  Figure 
14.  A bar  target  in  the  front  focal  plane  of  holographic  lens  was 
illuminated  with  light  made  spatially  incoherent  by  a rotating  diffuser 
(used  in  order  to  eliminate  speckle  in  the  image) . An  aperture  of  di- 
ameter 44  mm  was  placed  in  the  back  focal  plane  of  and  in  the  front 
focal  plane  of  the  reimaging  lens  L^,  which  was  designed  as  a colli- 
mating lens.  The  orientations  of  and  were  such  that  diffraction- 
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limited  imaging  was  expected  for  the  point  at  the  center  of  the  image. 

The  ratio  of  the  focal  lengths  of  and  being  300  mm/500  mm,  gave 
a net  magnification  of  0,6.  Figure  15(a)  shows  the  resulting  image. 

Figure  15(b)  shows  a magnified  portion  of  the  center  of  the  image, 
where  the  finest  resolved  bar  separation  at  the  input  target  is  seen 
to  be  about  7 pm,  which  is  roughly  equal  to  the  diffraction  limit  of 
1.22  AF/44  mm  = 7.1  ym.  Figure  15(c)  shows  a magnified  portion  of  one 
edge  of  the  image,  where  the  finest  resolved  bar  separation  is  about 
8 pm  vertically  and  11  pm  horizontally.  Assuming  a perfect  Fourier 

i 


transforming  lens  for  L^,  the  HOAD  program  predicted  a geometrical 
spot  of  rms  diameter  28  pm  for  that  location  in  the  image,  consider- 
ably worse  than  that  obtained  experimentally.  This  resulted  from  the 
fact  that  the  collimating  lens  L^,  a doublet,  does  not  perform  a per- 
fect Fourier  transform,  and  its  aberrations  partially  compensate  for 
the  aberrations  of  the  holographic  lens. 


VI.  CONCLUSIONS 

Analysis  was  performed  and  experiments  verified  that  holographic 
optics  can  perform  well  as  the  Fourier  transforming  elements  of  a 
matched-filter  optical  processor.  The  aberrations  of  a holographic 
before-the-lens  processor  were  found  to  be  comparable  to  those  of  an 
af ter-the-lens  processor.  Holographic  optics  of  a simple  design  are 
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capable  of  processing  imagery  of  large  space-bandwidth  product  if  long 
enough  focal  lengths  are  allowed.  The  standard  deviation  phase  error 
criterion  was  found  to  be  a useful  guide  for  predicting  the  performance 
of  the  lenses  in  a matched-filter  optical  processor. 

This  work  was  supported  by  the  Air  Force  Office  of  Scientific 
Research  under  Contract  No.  FA4620-76-C-0047 . 
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Figure  2.  Computer  simulation  of  matched  filter 
output  with  aberrations.  (a)  Test  ob- 
ject, (b)-(f)  central  portions  of  aber- 
rated autocorrelation  outputs:  (b)  no 
phase  error  (uppermost  curve)  and  O = 

A/16  nth  order  phase  error,  n = 1 to  4; 

(c)  = A/8;  (d)  = A/4;  (e)  20%  cen- 

tral stop,  no  phase  error;  (f)  20%  central 
stop,  a.  = A/4. 


Figure  5.  Translation  of  the  input  resulting  in  o,  = A/8, 
vs.  focal  length  (F) , for  various  widths  (Dg)  of  th 
image  of  interest.  A 10°  offset  and  40  cyc/mm  spa 
tial  frequency  are  assumed,  making  10  mm  equivalent 
to  800  picture  elements. 


Figure  6.  Performance  of  holographic  lens  for  Fourier  spectrum  analysis, 
assuming  a X/8  at  40  cyc/mm. 
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Figure  11.  Scan  through  peak  of  output,  = 4 mm. 
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Figure  15.  Output  image.  (b)  Magnified  central  area; 
(c)  Magnified  right  side  area. 
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